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Abstract 

Qh We derive formulae for neutrino masses and mixing angles in a type I seesaw framework with an underlying 

A4, flavor symmetry. In particular, the Majorana neutrino mass matrix includes contributions from an A^ 
triplet, 1, 1', and 1" flavon fields. Using these formulae, we constrain the general A^ parameter space 
^ using the updated global fits on neutrino mixing angles and mass squared differences, including results from 

^ the Daya Bay and RENO experiments, and we find predictive relations among the mixing parameters for 

00 certain choices of the triplet vacuum expectation value. In the normal hierarchy case, sizable deviation from 

\Q maximal atmospheric mixing is predicted, and such deviation is strongly correlated with the value of ^13 in 

the range of ~ (8— 10)°. On the other hand, such deviation is negligible and insensitive to ^13 in the inverted 
mass hierarchy case. We also show expectations for the Dirac CP phase resulting from the parameter scan. 
Future refined measurements of neutrino mixing angles will test these predicted correlations and potentially 

> 

show evidence for particular triplet vev patterns. 

>< 
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I. INTRODUCTION 



With the independent discoveries of a nonzero ^13 from the Daya Bay [T] and RENO |2] cohabo- 
rations, and the supporting hints from the T2K |3J, MINOS [4J, and Double Chooz [S] experiments, 
we now possess the first complete experimental picture of the Pontecorvo-Maki-Nakagawa-Sakata 
(PMNS) mixing matrix. Following a recent global analysis of neutrino oscillation parameters from 
Ref. [6j (see also [7H5] and [lUHT^ ). we can summarize the experimental status to date as 

sin2 012 = 0.3201°:°}^ , (1) 
. 2 . / 0-49tn^ (Normal) 

sm 6*23 = S ^nn, ' (2) 

[ 0.53l[!;[!^ (Inverted) 

/ 0.026t[J;0°3 (Normal) 
sm 6'i3 = < , (3) 

[ 0.027^[]:[!j]^ (Inverted) 

Am|i = 7.62 ± 0.19 x 10"^ eV^ , (4) 
2 f 2.53;0-?o X 10"' (Normal) 

^"T-Sl = \ ■ (5) 

[ -(2.40lg:J?) X 10"' eV^ (Inverted) 
Before the ^13 7^ discovery, the PMNS matrix was consistent with the tribimaximal (TBM) 
mixing pattern [13], which can be written as 

\/6 ^/3 
1 1 1 

y/Q yfi y/2 
1 1 1 

\ \/q \/3 \/2 y 

where we have adopted the same phase convention as the PDG [H] for placement of the minus 
signs. The TBM mixing pattern gives the solar mixing angle corresponding to sin^ ^12 = 1/3, 
the atmospheric mixing angle sin^ ^23 = 1/2, and reactor mixing angle sin^ ^13 = 0. Neutrino 
mass matrices that give rise to TBM mixing have distinct invariants that can be traced to discrete 
symmetries such as the Klein Z2 x Z2 group or the symmetry group 6*4 |15H18j (see also, [19]). On 
the other hand, by introducing dynamical (fiavon) fields, the TBM mixing pattern can arise from a 
smaller underlying finite group, the tetrahedral group, A4. To be compatible with grand unification, 
a successful generation of appropriate lepton and quark masses and mixing angles can be achieved 
by considering T', the double covering group of A^. In an SU{5) grand unified model [20], the T' 
group also affords a novel origin of CP violation from complex Clebsch-Gordan coefficients [21] . 
while in a Randall-Sundrum model |22il23j . the T' flavor symmetry is simultaneously used to forbid 



Utbm 



(6) 



tree- level flavor-changing neutral currents (FCNCs). For reviews of the status of A/j^ and ^4 models, 
the tribimaximal and bimaximal paradigm, see |24H26j . 

The literature on the interrelation between TBM neutrino mixing matrices and finite group 
flavor symmetries is vast. With discovery of a nonzero ^13 from Daya Bay and RENO, however, 
the TBM prediction for a vanishing reactor angle is ruled out, calling into question the entire TBM 
paradigm. We emphasize, however, that the underlying flavor symmetries that naturally give rise 
to TBM mixing are nevertheless viable options for explaining the updated PMNS mixing pattern. 
In particular, we demonstrate that A4, when we include flavons that are not included in the usual 
TBM analysis, can readily accommodate a large value of ^13 and retain predictivity for 5, the 
PMNS CP-violating Dirac phase. 

Work before the Daya Bay and RENO results that focused on finite group symmetries and 
generating nonzero ^13 includes |27j . which studied corrections to TBM from higher dimensional 
operators and [28], which discussed renormalization group (RG) equations in see-saw models. More 
recent work looking to use higher dimensional operators to generate violations of the TBM scheme 
or neutrino phenomenology include |29l l30] . The authors of [31] conclude higher dimension oper- 
ators and RG effects are equally important for leptogenesis. In [32], it was found that the size of 
corrections to neutrino mixing sum rules coming from renormalization group running are small. 
Correspondingly, in [33], it was shown that a large ^13 value cannot be generated through running 
if ^13 starts at 0. In [34J, NLO and NNLO expressions were given for neutrino mixing angles in 
hierarchal mass scenarios with sequential dominance. 

Nevertheless, we note that it is possible, with the so-called Hilbert basis method [35j, to con- 
struct supersymmetric models where these higher dimensional operators in the holomorphic super- 
potential vanish. On the other hand, certain fiavon-induced corrections in the non-holomorphic 
Kahler potential |36j cannot be forbidden by any conventional symmetries. Hence, even in the 
original models where ^13 vanishes at leading order, once these Kahler corrections are properly 
taken into account, a sizable 613 can be attained that is compatible with the current experimental 
value [37] . 

Much of the recent literature has focused solely on accommodating non-zero ^13 with low-energy 
perturbations to the TBM matrix and ignore or only briefly discuss the possible underlying UV 
physics responsible for these deviations. In this vein, one popular parametrization was introduced 
in [3H], which introduced deviations from TBM values using s, a, and r for the solar, atmospheric, 
and reactor mixing angles, respectively. Further work along this line, such as parameter scans in 
this space of deviations, include [391417] and also [IHHST] . which used a similar approach but an 
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equivalent set of parameters. One drawback of such approaches is that complete UV flavor models 
generally predict additional flavor violating effects, which are not captured in these low-energy 
deviation studies Separately, TBM and BM mixing scenarios were studied for the case of 

nearly degenerate neutrino masses where loop corrections provide large effects in mixing angles |53) . 
Loop corrections to the Type I seesaw Majorana neutrino mass matrix leading to nonzero ^13 were 
also studied in [54J . Studies of generating nonzero ^13 from breaking fi-r symmetry include |55H57] . 

In this work, we adopt a top-down approach where we will constrain the UV parameter space 
by the low energy neutrino observables. While other A4 studies may have looked at the additional 
effects from 1' and 1" flavons contributing to the Majorana neutrino mass matrix, including |58H64] . 
they have generally considered the pattern of the vacuum expectation value (VEV) of the triplet 
flavon to be oc (1, 1, 1)^. Some earlier work has separately considered other forms of the triplet 
flavon, such as [65]. In addition, recent literature has focused on small perturbations from the 
(1,1,1)"^ triplet flavon structure as a mechanism for generating a nonzero ^13 [58 l I66H68] . where 
this small perturbation may arise from a vacuum misalignment correction. 

In contrast, our work considers the full parameter space of A4 flavons contributing to the 
Majorana mass matrix in the Type I seesaw. Thus we consider triplet flavon vevs that are markedly 
different from the usual TBM (1, 1, 1)"^ form simultaneously with the presence of 1' and 1" flavons. 
By looking at a completely general admixture of possibly flavon vevs, we can definitively constrain 
the entire A4 parameter space in this Type I seesaw model. We remark that for the charged 
lepton masses and the Dirac neutrino masses, we introduce the minimal field content to generate 
their mass matrices and only introduce the full flavon field content for the RH neutrino masses. 
Moreover, even though we consider a larger A4 parameter space than the earlier literature, we still 
retain predictivity, especially when the triplet flavon vev pattern preserves a subgroup of A4. 

Here we concentrate on the group A4, which is the smallest group that contains a triplet rep- 
resentation. Other groups that have been utilized include 0(2) and SO{2) symmetry [69j, A(3n^) 
and A(6n^) |70l[7l], Agg [72], Qg with three sterile neutrinos [73], and product groups of modu- 
lar Zn finite groups in various permutations [73ti83] . Other studies have focused on permutation 
symmetry 5*3 and 5*4 models |84H9Uj . Other mixing scenarios beyond TBM and related patterns 
include democratic mixing, which has been studied in \4:2 \ \53 \ \9T\. tetra-maximal mixing [92j, and 
anarchic mixing |93H96j . 

The paper is organized as follows. In Sec. |llj we briefly review the A4 flnite group symmetry 
and our type-I seesaw model implementation. In Sec. |III[ we present the results of our parameter 



scan of the A4 flavor vev space. We conclude in Sec. IV, An intermediate step of our calculation 
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3 3 1 1" 1' 



3 3 1 1' 1 



TABLE I: The A4 charge assignments of the Standard Model fields and A4 flavons. 



is presented in Appendix |A] 
II. THE Ai MODEL 

We construct a Type I seesaw model based on an flavor symmetry. We include three right- 
handed neutrinos Ni, which are Standard Model gauge singlets. These neutrinos transform as a 
triplet 3 under A^. We also assign the lepton SU{2) doublet L ~ 3, charged lepton SU{2) singlets 
eR ~ 1, ~ 1", Tn ~ 1'. These representations are summarized in Table |lj 

The Lagrangian for the leptons is 



where (/>£;~3, i?i>Ar~3, ry~l, and -0 ~ 1" are scalar fields which acquire vevs and break 

the ^4 symmetry at the scale A, and the couplings cn, Cj^, c^, and are real. We will not specify 
here the scalar potential to give the (j)E, 4>Ni Vi Xi aiid fields vevs, but instead we leave a study 
of the potential construction and vacuum alignment questions for future work. The Lagrangian 
includes the familiar Dirac masses for the charged leptons, new Dirac masses for the neutrinos, and 
a general A4-invariant Majorana mass matrix which includes all possible ^4 contractions. Thus, 
in contrast with early, pre-Daya Bay and RENO models that generated exact TBM mixing by 
including only the ~ (1,1,1)^ and r/ flavons, we include the x and ij) flavons and allow (/)jv vev 
to be less constrained: subcategories of our approach have also been considered previously |58H64] . 

The explicit form of the resulting mass matrices for the charged leptons and neutrinos is easily 
obtained from the invariants, which are reviewed, for example, in Ref. [[24j. We assume (j)E 
acquires a vev {(pE) = ^(1)0,0)"^, and thus after electroweak symmetry breaking whereby the 
Higgs acquires a vev v^, the charged lepton mass matrix is 





/ Ae \ 



LMlc = Lvh 



Xa 



e 



(8) 



\0 Xr J 
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where e = (e/j, fipt, tr)- The corresponding Dirac mass matrix for the neutrinos is simply governed 
by the A4 contraction of L and A^, giving 

/ Xnvh ^ 



XnHLN = LMdcN = L 



\ 



N . 



(9) 



Xnvh 
Xnvh 

For the Majorana mass matrix, we wih be more general and allow (/>jv to obtain a general 
vev pattern, cn{4>n) = ^{(pa, (pb, (pc)- In addition, we let c^(?7) = Ar/, and similarly for x and V) 
such that the vev parameters r], and ip are dimensionless and have absorbed their respective 
Lagrangian couplings. We have the Majorana mass matrix 

-14)0 + ip -|4 + x^ 

V-g'/'fc + X -l4>a + ri l^c + V' J 
Since is a gauge singlet, there is no connection between the Higgs vev and the mass scale 
A^/j. In particular, if we set A/j/j ~ 0{Aqux) ^ W^^ GeV, we exercise the seesaw mechanism 
to generate light neutrino masses. Moreover, we assume the A4 breaking scale A ~ O.lA/j/j. The 
block matrix for the neutrinos in the {ul, N)"^ basis is 



MnN^N = ArrN'^N 



1^ +V' 106 + x 



(10) 



Ml^ 



Mdc Mn 

which generates the effective neutrino mass matrix 

M^, eff = MdcM^^MI^ 



(11) 



(12) 



after the right-handed neutrinos have been integrated out. 

Now, in exact analogy with the Cabibbo-Kobayashi-Maskawa (CKM) matrix, the PMNS matrix 
arises when we express the weak charged current interactions in the lepton mass basis. The PMNS 
matrix is 



where 



M' 



VpMNS = UlUI , 



(13) 



(14) 



but from Eq. ([s]), Ul = I3, and thus the PMNS matrix is identified with the neutrino diagonaliza- 
tion matrix ut. The three parameters Ae, A^, and A,-, which are each rescaled by {(j)E)/A, are in 
one-to-one correspondence with the three charged lepton masses and are hence fixed. 



Having identified Vpmns = ul, we adopt the standard parametrization of tlie PMNS matrix 
given by 



VpMNS 



/ 1 \ 

C23 S23 
\ -S23 C23 J 



( 



ci3 siae 
10 



-i<5 \ 



Y -5136*'' Cl3 / 



V 



Cl2 Sl2 
-Si2 Ci2 

1 



\ 



diag(l,e*«i,e*«2) , (15) 



wliere c 



cos ( 



sin^jj, 5 is the Dirac phase, ^1 and ^2 are the Majorana phases, with 



< < tt/2, 0<6, 6,6 < 27r. 

We now have a solvable system of equations relating the neutrino masses, mixing angles, and 



phases with the flavon vevs. From Eq. (|12|), Eq. (|14|), and using u}, = VpmnSi we find the relation 



Ml, cff = 



PMNS ' 



(16) 



where is given in Eq. (10). Solving for and inverting, we get 



Mn = Ml,VpMNs{Mf^^)-^Vl^,j,sMDc . 



(17) 



Entry by entry, we have a system of six equations that can be solved analytically, which gives exact 
relations between the ^4 parameter space and the physical masses and PMNS mixing angles. This 



is manifestly symmetric, which is most easily seen from Eq. (10). A similar result using the mass 



entries of the Majorana mass matrix instead of the triplet and one-dimensional flavon vevs was the 
starting point of [63J. Our approach, in contrast, directly shows the correlation between different 
UV triplet flavons and low energy neutrino observables. 

Naively, we have a six-dimensional UV parameter space, given by the three components of the 
triplet flavon and each of the one-dimensional flavons. The flavon breaking scale A and the 
Majorana mass scale are unobservable and can be absorbed into the six flavon vev components. 
Now, although the triplet's vevs can be independent degrees of freedom, certain breaking alignments 
in flavon space preserve the Z2 or Z3 subgroups of and thus reduce the number of UV parameters. 
We will thus categorize our results according to subgroup preserving and subgroup breaking triplet 
vev patterns, which are listed in Table |TIJ Subgroup preserving vev patterns thus effectively have 
a four-dimensional parameter space, while subgroup breaking vev patterns have a six-dimensional 
space. 

In general, there are nine physical parameters which give rise to eight physical predictions: 
three mixing angles, three masses, and one Dirac CP phase, as well as two Majorana phases that 
combine to dictate the rate of neutrinoless double (3 decay. We will not discuss Oz^2/3 further in 
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Subgroup 


vev alignment 


^2 


(1,0,0), (0,1,0), (0,0,1) 




(-1,1,1), (1,1,-1), (1,-1,1), (1,1,1) 


breaking 


(0,1,1), (1,0,1), (1,1,0), (0,1,-1), (2,1,1), (1,1,2), (1,2,1), (1,-2,1), (1,1,-2), (-2,1,1) 



TABLE II: Listed vev alignments for that preserve a or Z3 subgroup of A4. 



this paper, and instead focus on the three angles, three masses, and the Dirac CP phase. We can 
thus see both our subgroup preserving and subgroup breaking categories are predictive: the four 
parameters in the subgroup preserving category, for instance, are over-constrained by the existing 
neutrino measurements, and the existence of a nontrivial solution reflects the suitability of the A/^ 
finite group as a possible flavor symmetry of the lepton sector. In addition, for both categories, 
the yet-to-be-discovered Dirac CP phase is predicted from our parameter scan, which we detail in 
the next section. 



A. Breaking Bimaximality, Analytic Results 



We present an analytic understanding of breaking bimaximality (i.e. deviations of ^23 from 
45°) in our general Type I seesaw construction. First, we calculate the effective neutrino mass 
matrix, and then analyze to the extent a ^23 = 45° rotation diagonalizes this matrix. Starting 
with Eq. ([9]), we write 

/l 0^ 



Dc 



^NVh 



1 
1 



(18) 



and we can recognize that 



Mdc- Next, since Uu is the diagonalization matrix of 



M,y^ eff) as established in Eq. (14), 11^ is also the diagonalization matrix of ^^g, since 

(^diag)-i ^ ^^^^j^^ ^^f^t)-i = (C/t)-i(M,, ,srHUur' = Uu{M,, ,s)-^Ul . (19) 



Thus, the seesaw mass matrix in Eq. (12) can be inverted to give 



(20) 



Recall that because our charged lepton mass matrix in Eq. ([8| is diagonal, we have the identity 
VpMNS = ul, SO the diagonalization matrix of Mj^^ eff) and by extension, M~^^^ is the PMNS 
matrix. 
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Recall that the PMNS matrix is composed of three Jacobi rotation angles, ordered as a ^23 rota- 



tion, a rotation, and then a 9i2 rotation, and the phase rotations, as we can see from Eq. (15). 
Also, recall that a Jacobi rotation operates on a 2 x 2 block of the matrix, such as 




(21) 



where the rotation angle is defined as 



tan 2(9 



26 



(22) 



In particular, for fixed c and a, the sign of 6 determines the sign of 9. Moreover, we know that a 
rotation matrix for a single angle a + (3 can be decomposed into first rotating by a and then by /3. 
(This, of course, does not commute with rotations about other axes.) Hence, we can rotate M~ 
by ^23 = 45° and then understand deviations from bimaximality by testing the remaining presence 
of off-diagonal entries in the (3, 2) and (2, 3) entries of the effective neutrino mass matrix. 



From Eq. (|9|) and Eq. ([10|), the RHS of Eq. ^ is 

/ 2 



m: 



off 



X2 ,,2 



3^a + ri 3s 



(23) 



3v^^ 



4>h + 



which as a trivial check, is still symmetric. After we perform the ^23 = 45° bimaximal rotation, we 
have 

V ^('^^ + '^^) + V2(^ + '^) U-^b + <Pc)-lx+li^ l{-<l>a + (pb + (pc)+V+lx + l'4^ J 

(24) 

So, if |(— (/>fe + (pc) — |x + ^"0 7^ Oi then the ^23 = 45° bimaximal rotation was insufficient to 
eliminate the (3, 2) and (2, 3) entries and an additional ^23 rotation is needed. (In addition, if the 
(2, 2) and (3, 3) entries are identical, then 6*23 = 45° is guaranteed.) Moreover, we see that the 
transformation (pj, o (pc and x ^ changes the sign of the (3, 2) and (2, 3) entries while leaving 
the (2, 2) and (3, 3) entries fixed. Hence, we can see that any viable solution characterized by 
a triplet vev of {4>a-,4>b,<Pc) and a particular set of r/, Xi and ip can be transmuted to a different 
solution characterized by {(pa, fpc (pb) and ry, ip, and x with an opposite sign of the deviation from 
023 = 45°. 
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ai 


c?2C?3 - (■S12S23 - e "'cizcaasis) (-C23S12 - e"^ci2Si3S23)] 


a2 




C?3'5l2 + (ci2S23 + e "^C23Sl2Sl3) (ci2C23 - e*''si2Si3S23)] 


as 


e"^'«^(s?3-C?3C23S23) 


bi 


(si2S23 - e"*''ci2C23Sl3) (si2S23 " e"'ci2C23Sl3) + C12C13 (C23S12 + e""'ci2Si3S23)] 


b2 




(ci2S23 + e *''C23S12S13) (ci2S23 + e"'c23Sl2Sl3) " C13S12 (C12C23 - 6 "^512513523)] 


b3 


e"''^'(c?3cl3-e"'*Ci3Si3S23) 


Cl 


C12C13 (-S12S23 + e~"'ci2C23Sl3) + (C23S12 + e"^'ci2Sl3S23) (C23S12 + e~*''ci2Si3S23)] 


C2 


e"'2?i [ci3Si2 (ci2S23 + e^*''c23Sl2Sl3) + (ci2C23 - e*'^'si2Si3.S23) (C12C23 - e~'*Si2Si3S23)] 


C3 


e"''«^(c?3si3-e"'^Ci3C23Sl3) 



TABLE III: The explicit components of a, b, c as function of mixing angles and the triplet (j). 



III. PARAMETER SCAN RESULTS FOR SUBGROUP PRESERVING AND SUBGROUP 
BREAKING TRIPLET FLAVON VEVS 



Since current neutrino experiments do not have sensitivity to individual neutrino masses, we 
constrain the low energy neutrino observables by fitting to three mixing angles, two Am?, and the 
cosmological constraint on the sum of absolute neutrino masses. Clearly, the remaining parameter 
space for 5 is the predictive relation from our parameter scan. (As stated before, we do not discuss 
the sensitivity to the Majorana phases from neutrino- less double /3 decay experiments.) We first 



take the system of equations in Eq. (17) and solve for the breaking vevs in terms of the neutrino 
observables. These solutions are presented in Appendix |A| Since we want to work from the top- 
down, however, we partially invert the system to solve for the neutrino masses and the A4 singlet 
vevs in terms of the A4 triplet vev and the neutrino mixing angles and the CP phase. We obtain 

t;2 . . 

rrii = -r^ia'Vc''€ijk)i(l)ab''c>€ijk + cpba^c^eiju + (pca'^V eijk)~^ , (25) 
Afiij 

where mi, 1112, m^ are the three light neutrino masses, i, j, k = 1, 2, 3, and ejjfc is the Levi-Civita 
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tensor. The three-component vectors a, b, c are found in Table 

This partial inversion is advantageous because, using the measured mixing angles, we can test 
individual A4 triplet vevs and determine the fit to the correct mass squared differences and the 
cosmological constraint. We reparametrize the mass squared difference constraints into a ratio 
A21/A31, where Aij = mf — m| and use the dimensionless ratio to constrain the dimensionless 



vevs in Eq. (25). The mass scale vfj K^ji is then fixed by matching either of the measured mass 
squared differences. 

In line with our top-down approach and intuition about the A4 finite group, we attempt to pre- 
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serve many of the symmetries present in the case of TBM mixing. Yet, simply relaxing sin^ ^13 = 
while simultaneously keeping the other TBM mixing angle relations did not lead to viable solutions 
without perturbing the triplet vev alignment. We therefore choose to relax the bimaximal relation 
sin^ 623 = \ and maintain the trimaximal relation sin^ 612 = . In this way, we can understand a 
larger region of Ai^ parameter space since the experimental bounds on 612 are tighter than those 
on 023- Given this fixed ^12, we then choose each of the Dirac and Majorana phases to be or vr 
for each triplet vev pattern. 

However, this process itself is not entirely trivial, as not all choices are independent. Up to field 



0,6 


= 0,6 


= 5 = 


vr 


ei 


= 7r,6 


= vr 


0,6 


= vr,6 


= (5 = 


vr 


6 


= vr,6 


= 


0,6 


= 0,6 


= vr f-^ (5 = 


vr 


6 


= 0,6 


= vr 


0,6 


= vr,6 


= vr (5 = 


vr 


6 


= 0,6 


= , 



re-phasing, there are four transformations about bimaximal mixing in which components of Eq. ( 25 ) 
transform antisymmetrically. Under the interchange of 622, = 45° + x o ^23 = 45° — x and any of 
the four, 

(26) 



5 

5- 
we find 

a-{b'Kc)^ -a-{hxc) (27) 
6 X c •(-)• c X 6 
a X 6 o a X c . 

It immediately follows that the existence of one solution implies there is a corresponding solution 
with different phases (and a possibly different triplet alignment) which is related under these 
transformations. What remains to be specified are the angles ^23 and ^13, which we scan over the 
2a 623 range of [12j and over the 2a 613 experimental bounds from Daya Bay in Ref. [Ij . We 
draw contours in this two-dimensional plane that satisfy the measured mass squared differences 
from [I2J and we constrain the ratio of Am^ within 2a uncertainties as dictated by the predicted 
neutrino hierarchy and set the individual neutrino mass scale using A/j/j. From the existence of 
matter effects in solar neutrino oscillation, we know that A21 > 0. So by enforcing this constrain 
upon the ratio there is a unique positive value and negative value for A21/A31 that correspond to 
a normal or inverted neutrino mass hierarchy respectively. We then require the contours to satisfy 
the Xll^i^il < 0.81 eV [97j. Because degenerate neutrino spectra lead to numerically unstable 

i 

results, we discard possible solutions arising from degenerate neutrino spectra. 
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FIG. 1: Contour plots of A21/A31 over the 2a regions of ^23 and ^13 for subgroup preserving triplet vev 
patterns that generate (left panel) a normal hierarchy and (right panel) an inverted hierarchy. For normal 
hierarchies, we find the (1, 0, 0) pattern stays close to the central angle of 6'23, while the the (1, —1, 1) and 
(1, 1. —1) patterns avoid the central value. The inverted hierarchy solutions, however, only deviate from the 
maximal angle of 6*23 by values of order 0(0.1°). 

A. Numerical Results 



These results are shown in the left panel of Fig. [T] for subgroup preserving triplet vev patterns 
from Table [n| that generate a normal hierarchy, while the right panel shows the subgroup preserving 
patterns that generate an inverte(i hierarchy. These figures clearly show the correlations between 
the deviation from the bimaximal ^23 and nonzero ^13. Patterns that have such solutions are now 
uniquely determined and predict a very specific combination of masses and the CP phase. A 
notably absent alignment is the (1,1,1) breaking pattern, often associated with TBM mixing. We 
can understand this by considering Eq. (20) and asserting (pa = (f'b = <Pc- We find 



12, 1713 



,d23){m2 - m-i) = 



(28) 



and so unless the function / is zero, the masses m2 and m^, are required to be degenerate and 
this vev is not a viable phenomenological alignment. Under exact TBM mixing, / is zero and the 
masses are free. 

We also show the corresponding parameter scans for subgroup breaking triplet vev patterns 
that give normal hierarchies in the left panel of Fig. [2] and inverted hierarchies in the right panel. 
Although the parameter space is less predictive, we can nevertheless see that a normal hierarchy 
requires significant breaking of bimaximality, while an inverted hierarchy still allows for a bimaximal 
6*23 for some vev patterns. 
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FIG. 2: Contour plots of A21/A31 over the 2a regions of 6*23 and 6'i3 for subgroup breaking triplet vev 
patterns that generate (left panel) a normal hierarchy and (right panel) an inverted hierarchy. The nor- 
mal hierarchy solutions characteristically avoid the central value of ^23- Like the subgroup preserving vev 
patterns, here the inverted hierarchy solutions feature little deviations from the central value of 623 

For each successful triplet vev pattern and phase choice, we hst the 623 value and neutrino 
mass spectrum corresponding to the central value for 613 from the Daya Bay collaboration [l] 
in Table HVl 

The vev alignments (1, 0, 0) and the permutations of ( — 1, 1, 1) comprising the first and second 
sets of solutions shown in Table [IV| are vevs preserving the Z3 and Z2 subgroups, respectively. The 
remaining vev alignments break ^4 completely. Each vev alignment can produce some deviation 
from 623 = '7r/4, and most notably, non-Zs preserving alignments can produce significant breaking 
of bimaximality where TBM cannot be approximate. 

Generally, we find that inverted hierarchies feature deviations from 623 from the TBM case 
of order 0(0.1°), while normal hierarchies favor deviations an order of magnitude larger. For a 
given choice of Dirac CP phase, we find that most breaking patterns are relatively insensitive to 
changes in ^13 with the current experimental bounds. We also highlight the fact that many simple 
A4 triplet vev possibilities are excluded by the current experimental data. In particular, arbitrary 
triplet vev patterns in general will not generate the small hierarchy in experimental mass squared 
differences, and those that do are highly predictive about the extent of bimaximal breaking. 

We can semi-analytically see why the normal hierarchy solutions deviate more strongly from 
bimaximality than the inverted hierarchy solutions. In general, for fixed 9i2 and ^13, we can 
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consider a Taylor expansion of the neutrino masses around z = 623 — 45°, 

mi = xi + yiz + . . . , (29) 
ni2 = X2 + y2z + ■■■ , 
1713 = xs + ysz + . . . , 

where xi^2,3 and 2/1,2,3 are the zeroth and first order coefficients. Using this expansion, we have 

A21 =ml-ml^xl-xj + 2(x2y2 - xiyi)z + 0{z^) . (30) 

A complete expression of this mass squared difference for the central values of 9i2 and ^13 and 
arbitrary triplet vev is not useful, but if we consider the special case of 4>b = (j)c, we generally find 
xi = X2 and yi = -7/2, giving 

A21 ~ 4x22/22 + 0{z^) . (31) 



Furthermore, if we consider the behavior of the expansion coefficients xi and 2/1 for normal hierarchy 
solutions versus related in\ 
the A21 constraint, we see 



solutions versus related inverted hierarchy solutions, we find x\^^ ~ x\^^ and 4:y[^^ ~ Vi^^'- from 



z(^) « 4z(^) , (32) 

which implies that the resulting deviation from bimaximality is much larger for normal hierarchies 
compared to inverted hierarchies. When exact expressions are used, the difference can be as big as 
an order of magnitude, as evident in the figures. 

B. Dirac CP Phase Predictions 

We now examine the predicted Dirac CP phase for various triplet vev patterns. Instead of fixing 
6 and scanning over ^23 and ^13 as before, we now set ^23 to the value preferred by the central 



613 value of Daya Bay as shown in Table IV and scan this ^23 slice of the 6 vs. ^13 plane. 



Contours that satisfy the correct mass squared differences are highlighted and shown in Fig. [3] for a 
few illustrative choices of vev patterns. Thus, for ^13 within the 2a range of [Ij , the favored range 
of 5 can be broad, as for the (1,1,-2) vev, or fairly narrow, as for the (1,1,-1), (—1,1,1) and 
(0, 1, 1) vevs. This figure shows that a measurement of 6 and further refinement in narrowing the 
9i3 uncertainties can exclude or significantly favor a particular set of A4 vev patterns. In addition, 
we see that shifts in the central value of ^13 will serve to disfavor particular vev patterns as well as 
better accommodate other vev patterns. Certainly, more data is needed to test these possibilities 
and the paradigm. 
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IV. CONCLUSIONS 



In light of the results from the Daya Bay and RENO collaborations, we have developed a 
framework for understanding the constraints on the A4 parameter space from low energy neutrino 
observables. For our parameter scan, we have assumed a Type I seesaw model with a minimal 
A4 flavor structure governing the charged lepton masses and the Dirac neutrino masses, but the 
Majorana mass matrix has contributions from each of the possible 3, 1, 1', and 1" flavons, where 
the triplet vev pattern is initially unconstrained. We categorize the triplet vev according to Z3 
or Z2 subgroup preserving patterns, which enhances the model predictivity: for these subgroup 
preserving patterns, four A4 parameters are used to predict seven neutrino observables. Regarding 
the chosen vev patterns, we leave the question of scalar potentials or vacuum alignment for future 
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(1,-1,1) 
(1,1,-1) 
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(^,7r, 0,-4.8°) 
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I 
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(^,71,0,-0.2°) 
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N 
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TABLE IV: The collection of vev alignments considered in our parameter scan, resulting hierarchy, phases, 
deviation from bimaximality, and neutrino masses. The mass hierarchy column indicates whether the vev 
generates a normal (N) hierarchy or inverted (I) hierarchy. The last two columns indicate the required 
phases to obtain a valid mass hierarchy and the 623 angle and neutrino masses corresponding to the central 
9i3 value from Daya Bay [1]. 
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— (1,1-1) 
■ ■ (-1,1,1) 

■mill (0,1,1) 

T 8° 9° 10° 

FIG. 3: Contours for various A4 breaking patterns at their respective favored (?23 angle based on Table [IV| 
The (1, 1, —1) and (1, 1, —2) patterns show the ratio contour for a normal mass hierarchy, while the ( — 1. 1, 1) 
and (0, 1, 1) show the ratio for an inverted mass hierarchy. The (1, 1, —1), (—1, 1, 1), and (0, 1, 1) vevs are 
moderately more predictive of 5 than the (—1, 1, 1) vev. 

work. 

We have analyzed the parameter space in two distinct and intriguing slices. In the first case, 
we fix 9i2 to be trimaximal and scan over the resulting ^23 vs. ^13 plane for various choices of Dirac 
and Majorana phases. The results show that vev patterns giving a normal neutrino mass hierarchy 
have moderately large breaking of bimaximality for ^23) while inverted hierarchies generally retain 
bimaximality as a prediction. This indicates that a non-bimaximal ^23 measurement, such as the 
preliminary result sin^(2^23) = 0.94^qq5 from MINOS |98j is favorably correlated with a normal 
hierarchy in our framework. We also analyzed the predictions for the Dirac CP phase 6 for 
some illustrative choices of triplet vev pattern in the 5 vs. ^13 plane. This analysis emphasizes 
the point that a future measurement of (5, decreased uncertainty in ^13, and shifts in the central 
value of 013 can strongly favor or exclude particular triplet vev patterns, highlighting the fact that 
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future experimental results have significant power in discriminating possible ^4 flavon structures. 
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Appendix A: General VEV Solutions of Flavon Fields 



In Eq. (17), we have six equations that relate the A4 flavons comprising the Majorana mass 
matrix to the low energy neutrino masses, mixing angles, and phases. We can solve this system to 
express the triplet flavon components, {4>a-,4>b,(t'c)'^ , and the one-dimensional flavons, rj, x, and ip, 
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in terms of the physical neutrino observables. This solution set is 
1 
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(A6) 



where -F = ^ ^ . We invert Eq. (Al), Eq. (A2), and Eq. (A3) to solve for the neutrino masses 
^RR 



in terms of the triplet vev components {(f>a,4'b,4'c)'^- Then, we constrain the one-dimensional 
flavons and the neutrino masses by assuming a particular triplet vev pattern and scanning over 
the mixing angles and phases, applying the constraint on mass squared differences from and 
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the cosmological bound on the sum of absolute neutrino masses from |14] . 
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